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Using a combination of analytical techniques and Quan- 
tum Monte Carlo simulations we investigate the coupled spin 
ladder system LaCu02.s. At a critical ratio of the interlad- 
der to intraladder coupling (J'/J) c « 0.11 we find a quan- 
tum phase transition between a Neel ordered and a disor- 
dered state. At criticality the uniform susceptibility behaves 
as x{T) = olT 2 with a universal prefactor. At intermediate 
temperatures the system crosses over to a "decoupled ladders 
regime" with pseudo-gap type behavior, similar to uncoupled 
ladders. This can explain the gap-like experimental data for 
the magnetic susceptibility of LaCu02.5 despite the presence 
of long range Neel order. 

PACS numbers: 75.30Kz, 75.10. Jm, 75.40. Mg, 75.50.Ee 



The unusual normal state magnetic properties of doped 
high-T c cuprates have led to enhanced interest in zero 
temperature order-disorder transitions of quantum mag- 
nets. In particular, detailed predictions have been made 
about the behavior of a two-dimensional (2D) Heisenberg 
antiferromagnet by mapping it to the nonlinear sigma 
model 0]. They are in good agreement with experimen- 
tal measurements on La2Cu04. In addition to various 
mechanisms proposed for 2D spin systems, long range 
Neel order at T — can also be destroyed if a 3D an- 
tiferromagnet approaches the ID limit due to spatially 
anisotropic exchange. Then, quantum critical behavior 
and a disordered spin-liquid phase should be observed in 
three spatial dimensions. 

Recently, a suitable system for such type of behavior, 
LaCuC>2.5, has been synthesized The copper atoms 
in this compound form an array of coupled spin-1/2 two- 
chain ladders. Isolated spin ladders have a spin-liquid 
ground state and show signs of superconducting pairing 
with a d-wave order parameter upon doping ||. How- 
ever, a marked transition to a metallic phase takes place 
in LaCu02.5 under Sr doping, but no sign of supercon- 
ducting pairing was observed down to 5 K Q]. In con- 
trast superconductivity was recently found in the ladder 
compound Sro.4Ca13.6Cu24O41.84 which has weak and 
frustrated interladder couplings. This observation makes 
it quite important to study in detail the influence of the 
interladder coupling on the magnetic properties of the 
undoped insulating phase. 

First susceptibility measurements on LaCu02.5 were 
interpreted as showing a spin-gap in the excitation spec- 
trum |1. Subsequent NMR and /iSR studies indicated, in 
contrast, antifcrromagnetic ordering below T/v ~ 110 K 
HH ■ Normand and Rice suggested that the magnetic 



state could be close to a transition to spin-liquid phase. 
In this letter we expand on this idea and show that the 
apparently conflicting experimental results can be recon- 
ciled. 

The basic model for understanding these properties of 
LaCuC>2.5 is a spin-1/2 Heisenberg Hamiltonian for cou- 
pled ladders 
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which are shown schematically in Fig. [I]. We assume 
for simplicity equal rung and leg exchange constants J 
in each ladder and different exchange J' between lad- 
ders. Notice that the crystalline structure of LaCuC>2.5 
is more complicated, having four spins per unit cell [||. 
However, we may choose a simpler, topologically equiva- 
lent lattice structure having only two spins per unit cell. 
For J' w J the spin system is three-dimensional and has 
Neel order at low temperatures because the interladder 
coupling does not introduce frustration. Quantum fluc- 
tuations become more and more significant as one ap- 
proaches the quasi ID limit. Since the ID phase is a spin 
liquid with a finite gap, the magnetic order is destroyed 
at some finite J'. 

We examine the following points: (i) the critical ratio 
of (J'/J) c for the order-disorder transition, (ii) the low- 
T behavior of the uniform susceptibility x at the critical 
point, and (iii) x(T) in the whole temperature range and 
for arbitrary J'/ J. For this we employ a combination of 
analytical and numerical techniques. With the help of the 
renormalized spin- wave theory |g] and the bond-operator 
method Q we obtain lower and upper bounds for the 
transition point: 0.05 < (J'/J) c < 0.12. The quantum 
critical behavior of the uniform susceptibility for a 3D 
spin system has been predicted from scaling arguments 
by Chubukov et al. as X (T) = aT 2 . We calculate for 
the first time a universal factor in this law. Employing a 
Quantum Monte Carlo cluster algorithm (QMC) jl0| we 
then obtain a better estimate for the critical coupling: 
0.11 < (J'/J) c < 0.12. Next we calculate the temper- 
ature dependence of the uniform susceptibility x{T) for 
the whole temperature range and various coupling ratios, 
shown in Fig. |2j. Finally we show that the the suscepti- 
bility measurements of Hiroi and Takano can be fitted 
perfectly by the predicted form for a nearly critical or- 
dered system, thus resolving the apparent contradiction 
between the susceptibility and magnetic resonance mea- 
surements. 
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A natural approach to the Hamiltonian ([j]) from the 
side of strong interladder coupling J' ~ J is the renor- 
malized spin- wave theory of antiferromagnets Fol- 
lowing a slightly different procedure, we express the two 
spins per unit cell via two types of boson operators <ij 
and bi using the antiferromagnetic Dyson-Maleev trans- 
formation. Interaction terms with four bosons are then 
treated in the mean-field approximation by introducing 
boson averages: m = (afa^, Ai = (aA), A 2 = (aia t+z ), 
A3 = (aibi+ x ), which are determined by solving self- 
consistent equations. The corresponding spin-wave spec- 
trum consists of two branches 

cu k = ^4 2 -(# k ±|Ck|) 2 , (2) 
A = J{S-J7i+A 1 ) + 2J(S-m + A 2 ) + 2J'{S-J7i+A 3 ) 1 
By = 2J(S-m + A 2 )cosk z , 

C k = J(S-m+A 1 ) + 2J'(S-m+A 3 )(e lk * + e ik y) , 

each having zero-frequency mode at k z — or n. At 
the isotropic point J' = J, our calculations predict for 
S = 1/2 only a small reduction of the sublattice magne- 
tization: (S) = 0.40. Quantum fluctuations destroy the 
magnetic order for the critical coupling J' c 0.05 J. (The 
result by the linear spin- wave theory is an order of mag- 
nitude smaller.) From general arguments we expect that 
the renormalized spin-wave theory overestimates the sta- 
bility region of the ordered phase and, hence, J' c « 0.05 J 
presents a lower bound for the exact critical value. 

In the ordered phase the uniform magnetic suscepti- 
bility becomes anisotropic with two components paral- 
lel and perpendicular to the staggered moments. The 
parallel component x vanishes at T = 0. We calcu- 
late its low-temperature behavior in the framework of 
the present approach by using 

X I! = ^I>^>- (3) 

j 

In the limit T — > we find in agreement with Oguchi's 
results @]: x" = T 2 /Qc^c 1 ^, where C11 and c± are the two 
spin- wave velocities determined from (||) . The numerical 
coefficient in the square-law behavior of x(T) increases 
by a factor of 20 between J' = J and J' c . 

Describing correctly transverse oscillations in the or- 
dered phase, spin- wave theory fails, however, in the vicin- 
ity of J' c since at the critical point excitation spectrum 
has the same triplet degeneracy as in the disordered sin- 
glet phase for J' < J' c . To study the order-disorder tran- 
sition from the opposite side, we use the bond operators 
formalism This method describes a single spin lad- 
der fairly well for strong enough rung coupling [Q . It has 
also been applied to a 3D array of ladders in LaCuC>2.5 
at T = 0, but the result of Ref. JjJ is different from ours. 

The two spins (n = 1,2) belonging to the same ladder's 
rung with the lattice index i are expressed in terms of 
dimer states as 



(—\\ n j 
SZ ti = ^-{s+t^ + - -e^t+J^, (4) 

where Sj and t a i are singlet and triplet boson operators 
subject to the constraint sfsi + J2 a ^ai^a i = !• This re- 
lation is enforced by a chemical potential /1. Also, a site 
independent condensate of singlets (sj) = s is assumed. 
In the quadratic approximation we keep only the terms 
with two triplet operators. Diagonalizing the remaining 
Hamiltonian by the Bogoliubov transformation we ob- 
tain two self-consistent equations (<97i qU ad/<9/i) = and 
(&Hqaad/9s) — on the parameters \i and s. They can 
be reduced to a single equation on the new parameter 
rf = 2Js 2 /(J/4-/i): 

where = cos k z — J'/2 J(cos k x + cos k y ), magnon dis- 
persion is c^k = (>//4 — /iWl + dvy, and is a Bose 
factor. We first solve Eq. ([5|) at T = 0. The gap becomes 
zero for d = 1/(1 + J'/2J), Substituting this value into 
(||) we find that the critical coupling corresponding to 
vanishing gap and to the transition to the ordered phase 
is J' c = 0.121J. The mean-field theory should again over- 
estimate the stability region of the corresponding phase. 
Therefore, we conclude that the above value is an upper 
bound for the exact value of J' c , which lies between 0.05 
and 0.12. We will find below from QMC that the exact 
critical coupling is very close to the upper bound. The 
spectrum of low-lying excitations in the disordered phase 
near (0,0, 7r) has the form uiy = cm^/k| +p 2 k^ + m 2 , 
where C|| and c±_ = pc\\ are spin-wave velocities parallel 
and perpendicular to ladders, cy ~ 1.16J (at J' = J' c ), 
and p = y J'/2 J . The mass m and the gap A = c\\m 
behave like (J' c — J') 1 / 2 close to the critical point. 

The isotropic susceptibility in the spin singlet state can 
be calculated by Eq. (g), which after substitution of (^) 
takes the form 

*=^E( n k + «k), (6) 

k 

where summation is performed over one of the three 
magnon branches only. 

If the temperature is smaller than the gap, one can use 
the zero-temperature spectrum. In this quantum disor- 
dered regime the asymptotic behavior of the susceptibil- 
ity found from Eq. (JsJ) is 

A3/2 T l/2 

X(T) = * f 2 e~*' T , (7) 

which differs by its prefactor from the analogous results 
for magnetically disordered phases in ID and 2D 

At J' — J' c the mass m is generated by thermal fluctua- 
tions. It can be found from the self-consistency equation 
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at finite T. In contrast to the 2D case [12[, variation of 
the zero point fluctuation term in Eq. (|) becomes log- 
arithmically divergent on the upper limit, and is, there- 
fore, lattice dependent. Accordingly, m is a linear func- 
tion of T with logarithmically small prefactor computed 
by evaluating lattice sums: 



2n 2 



T 2 3 ln(0.7 J/T) 



(8) 



To calculate the universal behavior of the uniform sus- 
ceptibility in the quantum critical region A « T « J' 
we should neglect logarithmically small mass and substi- 
tute the gapless dispersion into Eq. (Q). As a result, the 
universal form for the susceptibility coincides with the 
result for x obtained in the spin-wave theory: 



X(T) 



T 2 



6c\\c 2 ± 



(9) 



Notice that nonuniversal corrections to the prefactor in 
the above expression have only logarithmic smallness. 

Analogous calculations for the specific heat predict 
C(T) = 2n 2 T 3 /5c\\c 2 L at the critical point. The tempera- 
ture dependence coincides again with the behavior in the 
ordered phase. However, the prefactor is multiplied by 
3/2 according to the different number of gapless modes in 
the two phases. Consequently, a crossover between these 
two regimes should exist for a "nearly critical" ordered 
spin system. 

Critical behavior can be also studied using a sigma 
model description of quantum antiferromagnets Q . Pre- 
dictions of that method have been compared with bond- 
operator results for a 2D magnet in Ref. . By analogy 
we argue that the limit N — > oo of the O(N) quantum 
nonlinear er-model in 3 + 1 dimensions should give the 
same universal factor as in Eq. (^|). This is quite natu- 
ral since both approaches use mean-field approximation. 
Calculation of leading 1/N corrections to the mean-field 
prediction remains an open question. 

Using QMC we can obtain a better estimate for the 
critical coupling. We have calculated the uniform sus- 
ceptibility x(T) for various couplings on lattices up to 
10 x 10 ladders of length 40 (8000 spins) and periodic 
boundary conditions at temperatures down to /3 J = 24. 
The results are shown in Fig. |[ We estimate the critical 
coupling by varying the coupling ratio and looking for the 
predicted T 2 behavior at criticality. Taking into account 
the shift of the critical point due to finite size effects |r| 
we estimate: 0.11 < (J'/J) c < 0.12, very close to the 
bond-operator estimate. 

Additionally we use self-consistent field boundary con- 
ditions H] to probe the occurrence of Neel order and to 
estimate Neel temperatures. We find Tjy sa 0.38(3) J at 
J' I J = 0.25, T N « 0.27(3) J at J'/ J = 0.15 and no indi- 
cation for order down to T = J/16 at J'/ J = 0.1. These 
results are consistent with the above estimates and show 



that the Neel temperature of about 110K (« J/10) ob- 
served in the experiments is realized very close to the 
critical point. 

Next we want to discuss x(^) f° r the whole temper- 
ature and coupling range and compare with the exper- 
imental measurements. For all couplings the Curie be- 
havior at high temperatures changes over into a broad 
maximum at temperatures of the order of J, caused by 
local spin singlet formation on the individual ladders, 
just as in uncoupled ladders [jll). The single ladder then 
shows a steep decrease with lowering the temperature, 
following an exponential decay x(T) ~ T _1 / 2 e _A / T pl| ] 
with a gap of about 0.5 J at low temperatures 

A weak coupling between the ladders does not destroy 
the spin gap. At high and intermediate temperatures 
we observe the same behavior and a steep exponential 
decrease with a pseudo gap similar to the gap of the 
single ladder. Only at temperatures of the order of J' 
a crossover to the 3D quantum disordered behavior Eq. 
(0) , an exponential decay with the actual gap, takes place. 

When A becomes smaller than J' upon approaching 
the T = transition point, the quantum critical region 
[jj] with its T 2 -law for the susceptibility appears between 
the quantum disordered and decoupled ladders regimes. 
Note, that existence of the 3D-type quantum critical be- 
havior is restricted to quite low temperatures T < J' . At 
T > J', when interladder coupling can be neglected, x(T) 
still shows a remarkable similarity to the single ladder. 

In the ordered phase close to criticality we find the 
same pseudo gap behavior, but the susceptibility goes 
to a small but nonzero value at zero temperature. The 
crossover occurs at temperatures of the order of the Neel 
temperature (compare J' / J = 0.15, 0.2 in Fig. |J). 

Let us now fit the susceptibility measurements on 
LaCuO 2.5. Hiroi and Takano have fitted them to an 
exponential form plus a Curie contribution due to impu- 
rity spins, and thus concluded a disordered ground state. 
But, as the magnetic resonance measurements indicate 
an ordered ground state the correct low-T behavior is 



X (T) = C/(T-e) + xo + aT 2 /j\ 



(10) 



where a w 0.33(3)e.m.u. mol -1 estimated from QMC, 
and xo is the sum of the temperature independent 
core susceptibility, and van Vleck susceptibilty and the 
small zero-temperature spin susceptibility. The fit is 
excellent, as shown in Fig. 3., with fitting parame- 
ters C = 1.8(1) x 10- 3 e.m.u.mol" 1 , 9 = -6.0(4)K, 
Xo = -6.2(2) x lO^e.m.u.mor 1 , and J = 1340±150K. 

We see that the uniform susceptibility measured by 
Hiroi and Takano || is indeed compatible with a gapless 
ordered ground state close to quantum criticality, as sug- 
gested by Normand and Rice 0. We remark that due 
to the dominance of quantum fluctuations in this nearly 
critical system no anomaly can be observed at the Neel 
temperature. 
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FIG. 1. Cross section of the lattice structure of the model. 
The ladders run perpendicular to the paper plane. Solid lines 
are the rungs of the ladders with a coupling J. Dashed lines 
are the inter ladder couplings J'. The dotted lines indicate 
the unit cell used. 

FIG. 2. Uniform susceptibility calculated by QMC for some 
representative ratios of the couplings. Error bars were omit- 
ted where the relative error was less than 1%. The inset is 
the same data in a double logarithmic plot. The dotted line 
is added as a guide to the eye, indicating the critical T 2 be- 
havior, clearer. 



Measurements of the total susceptibility suffer from the 
Curie contribution of impurity spins at low temperatures, 
which make the extraction of the asymptotic T — > be- 
havior difficult. Thus measurements which are not sensi- 
tive to impurities, such as NMR or /iSR are much better 
in distinguishing nearly critical ordered magnetic mate- 
rials from disordered ones. 
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FIG. 3. Fit of the susceptibility measurements by Hiroi and 
Takano |^] to Eq. (|lo|). The circles denote the measurements, 
the solid line the fit, and the full circles the measurements 
after subtraction of the Curie term. 
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Fig. 1, Troyer et al. 



Fig. 2, Troyer et al 
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Fig. 3, Troyer et al. 




